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Abstract. Mean-field models, while they can be cast into an extensive thermody- 
namic formalism, are inherently non additive. This is the basic feature which leads to 
ensemble inequivalence in these models. In this paper we study the global phase dia- 
gram of the infinite range Blume-Emery-Griffiths model both in the canonical and in 
the microcanonical ensembles. The microcanonical solution is obtained both by direct 
state counting and by the application of large deviation theory. The canonical phase 
diagram has first order and continuous transition lines separated by a tricritical point. 
We find that below the tricritical point, when the canonical transition is first order, 
the phase diagrams of the two ensembles disagree. In this region the microcanonical 
ensemble exhibits energy ranges with negative specific heat and temperature jumps at 
transition energies. These two features are discussed in a general context and the ap- 
propriate Maxwell constructions are introduced. Some preliminary extensions of these 
results to weakly decaying nonintegrable interactions are presented. 



1 Introduction 



Thermodynamics of systems with long range interactions is quite distinct from 
that of systems where the interactions are short ranged. For pairwise interactions, 
long range potentials decay at large distances as V{r) ~ 1/r" with < a < d in 
d dimensions. Such systems are not extensive , although models used to describe 
them can be made extensive by an appropriate volume dependent rescaling (the 
Kac prescription) |^ . However, even with such rescaling, these models remain 
inherently non- additive. Namely, when divided into two or more macroscopic 
subsystems, the total energy (or other extensive quantity) of the system, is not 
necessarily equal to the sum of the energies of the subsystems. Since additivity 
is an essential ingredient in the derivation of thermodynamics and statistical 
mechanics of systems with short range interactions, its violation may lead to 
interesting and unusual effects. 

To illustrate the lack of additivity in these systems consider an extremely sim- 
ple example of an Ising model with infinite range, mean-field like, interactions. 
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corresponding to a = 0. The Hamiltonian takes the form 




2 



where tr^ = ±1 is the spin 1/2 variable on site i and N is the number of spins 
in the system. This Kac type pre-factor is taken to ensure extensivity. Let us 
consider for simphcity a system with vanishing magnetization, M — X^i^i '^i ~ 
0. It is clear that the energy, E = H, oi the system satisfies E = 0. Now let us 
divide the system into two subsystems, / and II, each composed of N/2 sites, 
where all spins in subsystem / are up while those in subsystem II are down (see 
Fig. |l|). The energies of the two subsystems, Ej and Eji satisfy Ej — En = 
— JN/A. Since the sum of the two energies Ej + En is not equal to the total 
energy E the system is clearly not additive. 



I 


II 


+ 





Fig. 1. A two phase configuration of the infinite range interaction Ising model for which 
the energy is non-additive 



The lack of additivity may result in many unusual properties of systems with 
long range interactions, like inequivalence of various ensembles, negative specific 
heat in the microcanonical ensemble and possible temperature discontinuity at 
first order transitions. For example, the usual argument for the non-negativity of 
the specific heat, or equivalently, for the concavity of the entropy-energy curve 
s(e), makes use of the fact that systems with short range interactions are additive. 
Here s = S/N and e = E/N are the entropy and energy per particle, respectively. 
An entropy curve which is not concave in the energy interval ei < e < 62 (see 
Fig. 1^) is unstable for systems with short range interactions since entropy may be 
gained by phase separating the system into two subsystems with energy densities 
ei and 62 while keeping the total energy fixed. In particular the average energy 
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and entropy densities in the coexistence region is given by the weighted average 
of the corresponding densities of the two coexisting systems, and thus the correct 
entropy curve in that region is given by the common tangent line (see Fig. 
resulting in an overall concave curve. On the other hand for systems with long 
range interactions, additivity does not hold, and thus it is not always the case 
that entropy may be gained by phase separation. In such cases the non-concave 
entropy curve may in fact be the correct entropy of the system, resulting, for 
example, in negative specific heat. 




Fig. 2. A non-concave entropy curve, which for additive systems is made concave by 
considering two phase coexistence configurations in the region where the entropy of the 
homogeneous system is non-concave. This procedure is not applicable for non-additive 
systems like those with long range interactions. 



Self gravitating particles are perhaps the most extensively studied systems 
with long range forces (corresponding to a = 1 in d = 3 dimensions). The 
fact that the entropy of these systems is not necessarily a concave function of 
the energy has first been pointed out by Antonov Q. The thermodynamical 
consequences of this observation have been elaborated using simple models by 
Lynden-Bell |^ , Thirring Q and others (see also the papers by Chavanis Q and 
Padhmanaban ||^ in this book for a review). 

In fact, non-additivity is not limited to systems with long range interactions. 
Finite systems with short range interactions, such as atomic clusters, are non- 
additive as long as they are small enough so that surface energy may not be 
neglected as compared with the bulk energy . Negative specific heat in clusters 
of atoms has been discussed by R. M. Lynden-Bell js) and it has been recently 
observed experimentally, as discussed in this book in Refs. 
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In the present paper we consider a simple spin 1 mean-field model (a — 0) 
and study its phase diagram both within the canonical and the microcanonical 
ensembles. This model, known as the Blume-Emery-Grifhths (BEG) model p^ , 
may serve as a useful tool for studying thermodynamic effects which characterize 
systems with long range interactions. The canonical phase diagram of this model 
has been studied in the past, and it exhibits a line of phase transitions separating 
a magnetically ordered phase from a paramagnetic one. The line is composed of 
a first order transition and a second order transition segments separated by a 
tricritical point. The microcanonical phase diagram on the other hand exhibits 



some rather distinct features 1 11 . While the second order segment of the canon- 
ical phase diagram is present also in the microcanonical ensemble, the nature of 
the first order segment, its location and the location of the tricritical point are 
rather different. In this ensemble, a region with negative specific heat is found, 
as well as a first order transition at which the temperatures of the two coexisting 
phases are different. 

The paper is organized as follows: in Sec. || the BEG model is introduced 
and its canonical phase diagram is studied. The microcanonical phase diagram 
and the relationship between the two ensembles are studied in Sec. ^. Some fea- 
tures of the Maxwell constructions leading to temperature jumps are presented 
in Sec. ^. An alternative more general derivation of the expression for micro- 
canonical entropy is presented in Sec. ^. The results presented in this paper are 
generalized to models with slowly decaying interactions (with < a < rf) in 
Sec. 1^. Finally, concluding remarks are given in Sec. 



2 The Blume-Emery-GrifRths model in the canonical 
ensemble 

We consider a spin model with infinite range, mean-field like, interactions whose 
phase diagram can be analyzed analytically both within the canonical and the 
microcanonical ensembles. This study enables one to compare the two resulting 
diagrams and get a better understanding of the effect of the non-additivity on 
the thermodynamic behavior of the system. The model we consider is a simple 
version of the Blume-Emery-Griffiths (BEG) model [|lOj, known as the Blume- 
Capel model, with infinite range interactions. The model is defined on a lattice, 
where each lattice point i is occupied by a spin-1 variable Si — 0, ±1. The 
Hamiltonian is given by 

N T / ^ V 

i=i \j=i / 

where J > is a ferromagnetic coupling constant and A controls the energy 
difference between the magnetic {Si = ±1) and the non-magnetic {Si = 0) 
states. The canonical phase diagram of this model in the (T, A) plane has been 
studied in the past The free energy 

. 1 InZ 
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where Z is the partition function 

Z = ^exp-/3ff, (4) 

with /3 ~ X/ksT and ks the Boltzmann constant, can be exactly derived in the 
N ^ oo limit. One uses the Hubbard-Stratonovich transformation, which in this 
case amounts to the simple Gaussian identity 

— 1 dx exp{—bx^ + 2abx) , (5) 

I" J-oo 

with a — Si/N and b — (3JN/2. One then easily gets 



Z=\l^ I dx exp(- 7V/3/(/3, x)) (6) 



where 



x) = ^pJx^ - ln[l + e-^^(e'^^" + e-^'^^)]. (7) 

Using Laplace's method to perform the integral in formula @ in the iV — > oo 
limit, one finally gets 

/(/3)-min/(/3,.T) (8) 

At r = the model has a ferromagnetic phase for 2A/J < 1 and a non-magnetic 
phase otherwise separated by a first order phase transition. Indeed, the param- 
agnetic zero temperature state Si = 0,Vz is degenerate with the ferromagnetic 
state Si — l,Vi (or Si — — l,Vi) at 2 A/ J = 1 (this latter state being the ground 
state for 2Aj J < 1), a typical scenario for first order phase transitions. At 
A = —oo the model reduces to the mean-field Ising model and hence it has a 
second order phase transition at fc^T/J = 1. The (T, A) phase diagram displays 
a transition line separating the low temperature ferromagnetic phase from the 
high temperature paramagnetic phase. The transition line is found to be first 
order at high A values, while it is second order at low A. The second order line is 
found by locating the instability of the paramagnetic phase a; = 0, which means 
finding the condition for the x^ coefficient in formula (^ to be zero (note that 
x) is even in x). One gets for the second order term 

f3J = \eP^ + 1 . (9) 

The two segments of the transition line (high and low A) are separated by a 
tricritical point located at A/ J = ln(4)/3 ~ 0.4621, /3J = 3. This is obtained by 
requiring that both the x"^ coefficient of /(/3, x) in (0) and the x"^ coefficient 

1 2exp(-/3Z\) 

3 (l + 2exp(-/3Z\) ^ ' 
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vanish. The first order segment of the transition line is obtained numericaUy 
by equating the free energies of the ferromagnetic and the paramagnetic states. 
The behavior of the function /3/(/3, x) as /? varies is shown in Fig. ^ at a second 
order phase transition {A — 0) and in Fig. ^ at a first order phase transition 
{A = 0.476190). 



-1.1, 




Fig. 3. Pf{(3, x) as a function of (3J £ [1.4, 1.7] and x € [-1, 1] for A = Q. The second 
order transition is at (3, J = 3/2 where the minimum at a; = becomes a maximum and 
two side minima develop 



Microcanonical solution of the Blume-Emery-GrifRths 
model 



The derivation of the phase diagram of the BEG model (Q) in the microcanonical 
ensemble reduces to a simple counting problem, since all spins interact with equal 
strength, independently of their mutual distance (there is no space in the model). 
A given microscopic configuration is characterized by the numbers N+, Nq 
of up, down and zero spins, with iV+ + iV_ + N^) = N . The energy E of this 
configuration is only a function of iV+ , iV_ and iVo and is given by 



AQ 



J 



(11) 



where Q = X^iLi = ^+ + ^- (the quadrupole moment) and M = J^iLi '^^ ~ 
— N- (the magnetization) are the two order parameters. The number of mi- 
croscopic configurations fi compatible with the macroscopic occupation numbers 

N+,N^ and Nq is 

^2 = , . (12) 



N 



N+\N. 
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Fig. 4. Pf{(3,x) as a function of (3 J G [3.7,4.4] and x G [-1, 1] for A = 0.476190. The 
first order transition is at /3J = 4. Two side minima to x = first develop and then 
reach the same hight at the transition point. 



Using Stirling's approximation in the large N limit, the entropy, S = ksl'a f2, is 
given by 

S = -ksNiil - q) ln(l - q) + ^{q + m) ln{q + m) 

+ i(g-m)ln(g-m)-gln2] , (13) 

where q — Q/N and m = M/N are the quadrupole moment and the magnetiza- 
tion per site, respectively. Let e = E/ AN be the dimensionless energy per site, 
normalized by A. Equation ( pT| ) may be written as 

q = ( + Kr-n^ , (14) 

where K = J/2A. Using this relation, the entropy per site s = S/{kBN) can 
be expressed in terms of m and e. At fixed e, the value of m which maximizes 
the entropy corresponds to the equilibrium magnetization. The corresponding 
equilibrium entropy s(e) s(e, m) contains all the information about the 

thermodynamics of the system in the microcanonical ensemble. For instance, 
temperature can be obtained from the relation 

^ - - (15) 
keT - de ■ ^^^^ 

As usual in systems where the energy per particle is bounded from above, the 
model has both a positive and a negative temperature region: entropy is a one 
humped function of the energy. The interesting features take place in the pos- 
itive temperature range. In order to locate the continuous transition line, one 
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develops s{e,m) in powers of m, in analogy with what has been done above for 
the canonical free energy 

s ^ So + Am^ + Bm'^ + 0{m^) , (16) 

where 

So = s(e,m = 0) = -(1 - e) ln(l - e) - elne + eln2 , (17) 

and 

£ 1 



A = -K\n 



2(1 -e) 2e 



In the paramagnetic phase both A and B are negative, and the entropy is max- 
imized by m — 0. The continuous transition to the ferromagnetic phase takes 
place at A = for B < 0. In order to obtain the critical line in the (T, A) plane 
we first observe that temperature is calculable on the critical line (m = 0) using 
(© and (10) 

Requiring now that A = 0, one gets the following expression for the critical line 

2$K = + 1 , (20) 

where f3 = (3 A. Equivalently, this expression may be written as f3K = l/2e. The 
microcanonical critical line thus coincides with the critical line (||) obtained for 
the canonical ensemble. The tricritical point of the microcanonical ensemble is 
obtained at ^ = S = 0. Combining these equations with Eq. ( p^ ) one finds that 
at the tricritical point /? satisfies 

-i + i = 0. (21) 



8/32 e/3 4/3 12 

Equations (|o[ |l|) yield a tricritical point a,t K 1.0813, (3 ~ 1.3998. This has 
to be compared with the canonical tricritical point located at K — 3/ln(16) ~ 
1.0820, f3 = ln(4) ~ 1.3995. The two points, although very close to each other, 
do not coincide. The microcanonical critical line extends beyond the canonical 
one. In the region between the two tricritical points, the canonical ensemble 
yields a first order transition at a higher temperature, while in the microcanon- 
ical ensemble the transition is still continuous. It is in this region, as discussed 
below, that negative heat capacity appears. A schematic phase diagram near the 
canonical tricritical point (CTP) and the microcanonical one (MTP) is given in 
Fig. 1^. Beyond the microcanonical tricritical point the temperature has a jump 
at the transition energy in the microcanonical ensemble. The two lines emerging 
on the right side from the MTP correspond to the two limiting temperatures 
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Fig. 5. A schematic representation of the phase diagram, where we expand the region 
around the canonical (CTP) and the microcanonical (MTP) tricritical points. The sec- 
ond order line, common to both ensembles, is dotted, the first order canonical transition 
line is solid and the microcanonical transition lines are dashed (with the bold dashed 
line representing a continuous transition) 



which are reached when approaching the transition energy from below and from 
above (see Fig. ^d). The two microcanonical temperature lines and the canonical 
first-order transition line all merge on the T = line at 2A/ J = 1. 

To get a better undertanding of the microcanonical phase diagram and also 
to compare our results with those obtained for self-gravitating and for finite 
systems we consider the temperature-energy relation T(e). This curve has two 
branches: a high energy branch ( p^ corresponding to m = 0, and a low energy 
branch obtained from ( |l5| ) using the spontaneous magnetization ms{e). At the 
intersection point of the two branches the two entropies become equal. Their 
first derivatives at the crossing point can be different, resulting in a jump in 
the temperature, i.e. a microcanonical first order transition. When the transi- 
tion is continuous in the microcanonical ensemble, i.e. the first derivative of the 
entropy branches at the crossing point are equal, our model always displays, 
at variance with what happens for gravitational systems, a discontinuity in the 
second derivative of the entropy. This is due to the fact that here we have a true 
symmetry breaking transition. Fig. ^ displays the T(e) curve for increasing val- 
ues of A. For A/ J — ln(4)/3), corresponding to the canonical tricritical point, 
the lower branch of the curve has a zero slope at the intersection point (Fig. 
Thus, the specific heat of the ordered phase diverges at this point. This effect sig- 
nals the canonical tricritical point as it appears in the microcanonical ensemble. 
Increasing A to the region between the two tricritical points a negative specific 
heat in the microcanonical ensemble first arises {dT/de < 0), see Fig. |b. At 
the microcanonical tricritical point A the derivative dT/de of the lower branch 
diverges at the transition point, yielding a vanishing specific heat (Fig. ^). For 
larger values of Zi a jump in the temperature appears at the transition energy 
(Fig. ^). The lower temperature corresponds to the m = solution (^9|) and 
the upper one is given by exp(/3) = 2(1 — q*)/\/{q*)^ ~ (m*)^, where m*,q* 
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are the values of the order parameters of the ferromagnetic state at the transi- 
tion energy. The negative specific heat branch disappears at larger values of A, 
leaving just a temperature jump (see Fig. ||e). In the A/J 1/2 limit the low 
temperature branch, corresponding to q — m = 1 in the limit, shrinks to zero 
and the m = branch (19) occupies the full energy range (Fig. 




Fig. 6. The temperature-energy relation in the microcanonical ensemble for different 
values of A. The dotted horizontal line in some of the plots is the Maxwell construction 
in the canonical ensemble and identifies the canonical first order transition temperature 
at the point where two minima of the free energy coexist. The vertical line in panels d) 
and e) marks the presence of a temperature jump; for this energy two entropy maxima 
of the same height coexist 



We now present some general considerations concerning the intricate relation 
between canonical and microcanonical ensembles. The following discussion is 
quite general, and not specific to the BEG model. Let us first calculate the 
canonical partition sum using the density of states Q — e^'^^'^'"^\ where e = E/N 
is the energy density and m the order parameter (for the BEG model e = eA, 
m is the magnetization and q is obtained from e and m by relation (p^)). 

Z = ^exp(-/3i7(5,)) 

= J2 e^^f^^^^e-^'^^ , (22) 

[e,m] 

where we have replaced the sum over the configurations by the sum over the 
values of the energy per spin and the order parameter, which for finite N take 
discrete values. In the thermodynamic limit N oo, the free energy in the 
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canonical ensemble is thus given by 



/(/?)= min(e-^%^^) . (23) 



e.rn 



f3 



The point where the minimum is reached yields the mean energy and magnetiza- 
tion per spin in the canonical ensemble. As explained above, the microcanonical 
entropy is 

s(e) = max s(e, m) (24) 

m 

and the magnetization in the microcanonical ensemble is the point where the 
maximum is reached. 

All questions are now reduced to the solution of the two variational problems 
( p3| ) and (p^). The first order conditions are the same for the two problems 

dm 

where the second condition is indeed given in the microcanonical case by the 
definition of temperature. We will denote by e* {(3) , m* {f3) the solution of the 
variational problem (p5|). Using (p5|), it is straightforward to verify that 

'-^-e-Wl (26, 

meaning that the canonical mean energy coincides with the minimizing energy. 
Hence, the extrema in (|3|) and in ( p^ ) are the same. However this does not 
imply ensemble equivalence: we have to study the stability of these extrema. 

In order to discuss the stability of the canonical solution one has to determine 
the sign of the eigenvalues of the Hessian of the function to be minimized in (p3|). 
The Hessian is 

f ^mm ^7rte J (27) 

where, for example, Smm is the second derivative of s with respect to m. The 
extremum is a minimum if and only if the determinant and the trace of the 
Hessian are positive 

See Smm ^ (^^) 

> 0, (29) 



„2 



See^mm S^. 



which implies that See and Smm must be negative, and moreover See < —Sme/l^mm 
This has strong implications on the canonical specific heat, which must be posi- 
tive. Let us prove it using the variational approach, instead of the usual Thirring 
argument H, that uses the expression of the canonical partition sum. Indeed, 
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taking the derivatives of Eqs. ( p5| ) with respect to /?, after having substituted 
into them e*(/3), to*(/3), one gets 

de* dm* 

de* dm* , . 

where it is understood that all second derivatives are computed at e* {(3) , m* {j3) . 
Recalling now that the specific heat per particle at constant volume is, 

Cv de* ^,d^ ^32) 



NkB dT ^ di3 
one gets 

]Vfc7 ^ ^ .2 _ g , ' (33) 

which is always positive if the stability conditions (^) are satisfied. However, 
since the stability condition in the microcanonical ensemble only requires that 
Smm is negative, in order to have an entropy maximum, this indirectly proves 
that a canonically stable solution is also microcanonically stable. The converse 
is not true: one may well have an entropy maximum, Smm < 0, which is a free 
energy saddle point, with s^e > 0. This implies that the specific heat ( ^ ) can 
be negative . 

The above results are actually quite general, provided the canonical and micro- 
canonical solutions are expressed through variational problems of the type ( p3[ ) 
and (p^. The extrema, and thus the caloric curves T(e), are the same in the 
two ensembles, but the stability of the different branches is different (see also the 
important paper by Katz and the contribution by Chavanis to this book Q 
for a discussion of this point in connection to self-gravitating systems) . Another 
example of this behavior is discussed for the HMF model in this book p6[. 



4 Maxwell constructions 



In this Section we briefly comment on the Maxwell construction leading to tem- 
perature jumps in the microcanonical ensemble. This is quite a new feature, a 
sort of analogue of the latent heat phenomenon in the canonical ensemble. 

Let us begin with the ordinary Maxwell construction in the canonical ensem- 
ble. We consider the curve /3(e) as in Fig. 0a. It is the situation of Fig. ||c, if for 
simplicity we disregard the discontinuity in the temperature derivative, which is 
irrelevant for the present reasoning. The equal area Maxwell condition Ai = A2 
simply reads 

/" '[/3(e)-/3e]de = , (34) 

where (ic is the first order transition inverse temperature. Since (3 = ds/de, this 
implies 

s(e3)-s(ei)-/3c(e3-ei) -0. (35) 
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Using now the definition of free energy, f — e — s/ (3, we obtain the usual equal 
free energy condition at a first order phase transition, 



/(ei)-/(e3) 



(36) 



Near a discontinuous microcanonical transition, the /3(e) curve has the typical 
shape given in Fig. This is in fact the case of Fig. ||, if one continues the low 
energy branch of /3(e) above the microcanonical transition energy, drawing the 
metastable and the unstable lines as well. We invert this relation to obtain e(/3) 
in the vicinity of Cc. Then, the equal areas condition reads 



/33 



[e(/3) - ej d/3 = 



(37) 

(38) 
(39) 

This is the equal entropy condition at a microcanonical discontinuous transition. 
Eq. ( |37| ) is valid only for the simple configuration of Fig. |^; for more complex 
curves one has to evaluate the areas more carefully (see also the discussion by 
Chavanis Q). 



using e = d(/3/) /d/3, we get 

Mm - f3if{f3i) - ee(/33 -Pi) = 

which implies 

s(/33) = s(/3i) . 




/3" 

03 

Pi 







b) 










Ai) 





Cc e 
Fig. 7. Canonical (a) and microcanonical (b) Maxwell constructions. 
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5 The microcanonical solution by large deviation 
techniques 

In Section |^ we derived an expression for the logarithm of the density of states 
at fixed energy and magnetization, i.e. the entropy s(e,m), using combinato- 
rial techniques. We now present an alternative derivation of the entropy using 
Cramer's theorem This theorem, derived in the context of large deviation 
theory , states the following: consider a set of TV independent and identically 
distributed random variables in i?'^, (Xi)i=i...Ar. One would like to know the 
probability distribution of the average = 1/N^Xi. For this purpose, one 
defines, for X £ R'^, 

V'(A)=ln(e^-^), (40) 

where the average () is taken with respect to the common Probability distribution 
Function (PdF) of all Xi. Now let s{x)/kB, with x G R'^, be the inverse Legendre 
transform of ip 

4^ = - sup (X-x-^iX)) . (41) 

Then, under quite general conditions, Cramer's theorem states that the PdF of 
the average is 

PiMN ^x)^ e ~ . (42) 

This yields an expression for the logarithm of the density of states s{x)/kB- The 
great advantage of this method is, of course, its generality, since it applies also 
when combinatorial tools are not available. 

Let us give a brief heuristic argument for this result in the d ^ I case, 
A, X e R. The probability of having A/jv = a; is given by the volume in phase 
space compatible with Mn = x. Let dfi be the common measure of all Xi. One 
gets, 

P{Mn = 2;) = y dfi{Xi) . . . dfi{XN)S{MN - x) (43) 

= W-f rfAe-^^" / d^l{Xl)...d^^iXN)e^^^'+■■■+^''^ (44) 
Stti Jp J 

= -L / dXe-^'^[{e'^)f , (45) 

ZTTl Jp 

where the Dirac S function has been represented by a Laplace integral over a 
path r transverse to the real axis in the complex A plane. Evaluating the last 
integral using the saddle point method one is lead to look for the critical points 
of Ax - ^(A), with V^A) = ln[(e^-^)], which justifies Eq. (|l|). 

As an example, let us apply this method to the BEG modelQ In this case 
the Xi variables are bidimensional {d = 2), Xi = {Sf, Si), so that 

V'(A,p) = In (£^■^'+''■5') 



^ An application to the Ising model and to other systems can be found in 
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= ln(l + 2e^ coshp) - In 3 (46) 

To calculate s{q, m) we now have to solve equations dip/dX — q and dtp/dp = m 
for A and p. One gets 



where r is the ratio m/q. Substituting in (p|), we obtain the expected result 

s{q,m) q + m q + m q - m q - m n i o (An\ 

— ^— = ^ In — — In — {l-q) hi(l - g) - ln3 , (49) 

where the — In 3 constant is a normalization factor related to the fact that the 
total number of configurations is 3^. 



6 Slowly decaying interactions 

So far, we have restricted our study to infinite-range models. Space is in this 
case irrelevant, which simplifies very much the study. However, the physically 
interesting interactions (gravity for instance) are not infinite range, and one 
may wonder whether the properties of the previous Sections still hold in a more 
general case. To answer this question we discuss in this Section spin systems on 
d-dimensional lattices with pairwise interaction potential between two sites i and 
j decaying with the distance like ^/rf^ (0 < a < d), so that the interaction is 
long-range. An Ising model with short range algebraically decaying interactions 



(a > d) has been studied in the past by Dyson 1 15 . More recently, versions of the 
Hamiltonian Mean Field model , with slowly decaying long-range interactions 
among rotors (a < d) have been considered by several authors [ p7|p^ and the 
problem has been discussed in a more general framework in Ref. |19[. 

Our main result is that all peculiar features of the microcanonical ensemble 
of infinite range models extend to these slowly decaying interactions. We proceed 
by considering a continuum, or coarse grained version of magnetic systems. We 
then explain that in fact these coarse grained models can be derived from the 
microscopic models, like the Ising or the BEG models considered in this paper, 
using the large deviations tools introduced in Section 5 (see ||2^ for a more 
thorough discussion of this point). 

We begin by analyzing the phase diagram of the continuum version of the 
Ising model with algebraically decaying interactions (0 < a < c?) . We show that 
as in the mean field model {a = 0) this model does not exhibit a phase transition 
in the microcanonical ensemble. We then consider the continuum version of the 
BEG model. 

The Ising Hamiltonian functional with algebraically decaying interactions in 
d = 1 dimensions on the segment [0, 1] with periodic or free boundary conditions 
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(to be discussed below) takes the form 

J , , m{x)m{y) 



his.n, = -^J^ d.J^ dy^:^^, (50) 

where m{x) is the local magnetization, r{x, y) is the distance between x and y 
on the segment (the shortest distance for periodic boundary conditions), and 
J > is a ferromagnetic coupling constant . The entropy density corresponding 
to a magnetization profile m{x) of an Ising variable may be written as 

, /"^jrll+H, (1 + m) (1-m)^ (l-™)i 
SIs^ng = "^5 / dx ^- ^ In , + , hi , , (51) 



Jo A ^ ^ A 

Once the energy and entropy functionals are given, the microcanonical (canon- 
ical) solution is obtained by a maximization (minimization) of the entropy (free- 
energy) functional under the constraint of constant energy (temperature). We 
begin below with the discussion of these two variational problems for the Ising 
model 

s/smg (e) = max (s7s.mg[(m(a;)] I huing = e) , (52) 

m{x) 

in the microcanonical ensemble , and 



1 

i(x) \ p 



fistngiP) = min hising[m{x)] - -sising[m{x)\ ) (53) 



in the canonical ensemble . 

So far, although we have stated the problem for the Ising model, we could 
have followed the same path, whatever the lattice model and the boundary con- 
ditions are. We can reach a general conclusion at this point by writing down the 
extremality conditions of the two variational problems (^) and (|5^). The first 
order reads 

I sing p^hjg^ng (54:) 

in the microcanonical ensemble, where /3 is a Lagrange multiplier. In the canon- 
ical ensemble 

A/i, . 1 Ac, . 

. (55) 



6m{x) (3 Sm{x) 

Hence, the first order extremality conditions lead to the same equations for the 
two ensembles. However, this does not imply ensemble equivalence, since the 
stability of the solutions may differ in the two ensembles, as has been noted in 
Section |. 

Before proceeding, let us briefiy recall the results of Ref. |l^. These au- 
thors, discussing the Hamiltonian Mean Field model with slowly decaying inter- 
actions [ p^lp^ , remarked that, choosing appropriately the renormalization factor 
in the Hamiltonian, the thermodynamic behavior of the system turned out to 
be independent of a; i.e. the system behaves at < a < 1 (slowly decaying 
case) exactly as if a = (mean- field case). An explanation of this result was 
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provided by Campa et al. p8[ through a canonical analysis of long-range inter- 
acting systems on lattices with periodic boundary conditions, and independently 
by VoUmayr-Lee and Luijten in a more general context. Unfortunately, nei- 
ther of the two latter methods provide the microcanonical solution. The method 
we present here is, instead, fully general and allows to discuss in detail how to 
extend the results for a = to a > 0. 

Let us solve now the variational problem of Eqs.jjS^) and (|55| ) for the Ising 
model The magnetization profile that extremizes entropy and free energy 
satisfies the following self-consistency condition 



m(x) = tanh ( pj [' dy ] (56) 







In the canonical ensemble, (3 is fixed, whereas in the microcanonical ensemble 
it has to be tuned in order to get the required energy. For periodic boundary 
conditions, the integral on the r.h.s. of Eq. (Bq) does not depend on x, and 



Eq.(56) always has a homogeneous solution with 'm{x) independent of x. Eq. ( p6|) 
then reduces to the usual consistency equation of the mean-field Ising model, and 
the system behaves exactly as in the a = case, showing a phase-transition in 
the canonical ensemble (no transition is present in the microcanonical ensemble). 

On the contrary, for free boundary conditions, there is no non-zero constant 
solution, and the system develops a non trivial magnetization profile that can be 
only numerically computed from Eq . (^6|) . Moreover, in this case a phase tran- 
sition is present in the canonical ensemble, while the microcanonical ensemble 
is always in a magnetically ordered phase. The magnetization profile in the mi- 
crocanonical ensemble for free boundary conditions is shown in Fig. ^ for three 
difi^erent values of a. 




(1 i 

Pcisitiun ciii the lattitif 

Fig. 8. Magnetization profile for the Ising model with slowly decaying interactions 
in the microcanonical ensemble. The energy density is fixed at e = —0.1 while a 
varies: a = 0.2 (full line), a = 0.5 (dotted line), a = 0.8 (dashed line). The boundary 
conditions are free. 
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Let us now discuss the behavior of the BEG model for sfowly decaying inter- 
actions. The energy functional is 



= (57) 



where all is defined as before and, moreover, q{x) is the quadrupolar field in 
the continuum on the segment [0, 1] (obtained in an analogous way as for the 
magnetization m{x)). The entropy functional of the BEG model is obtained from 
expression ( p^ by replacing m, q by m{x), q{x). 

The BEG model is much richer, as has already been demonstrated for a = 0. 
The analysis of the microcanonical phase diagram for periodic boundary condi- 
tions follows the same path as that presented above for the Ising model. It is 
straightforward to show that the constant profiles m,{x) = m*, q(x) = g*, with 
m* and q* given by the solution of the a = model, are solutions of the Eqs. 
and ( |55| ) (where we replace the Ising subscript by the BEG one) also for a > 0. 
Moreover, detailed calculations show that these solutions are local maxima 
of the entropy and local minima of the free energy. Thus one concludes that 
for periodic systems the entropy and free energy are independent of a and the 
conclusions reached for a — Q are valid for < a < d. In particular, ensemble 
inequivalence can be extended to the case of long-range algebraically decaying 
interactions. 

Let us briefly sketch these calculations. To verify that uniform profiles are 
not destabilized by modulations, it is convenient to develop q{x) and m{x) in 
Fourier series: q{x) — '^f^qkUk(x), m{x) — '^f.rhkUkix), where k ranges from 
— oo to -l-oo and Uk(x) = coskx if fc > 0, Uk(x) = sinfca; if fc < . 

The Fourier representation, due to the periodicity, diagonalizes the Hamilto- 
nian. The energy may thus be written as 

e = qo-K ^kfhl . (58) 

fc— — oo 

Here Afc are the energy eigenvalues and K = 2" J/2Z\(1 — a) is chosen such that 
the largest eigenvalue is unity Aq = 1 (this amounts to choose the prefactor N 
in the Kac prescription, see below). 
The free energy functional is 



/[q, m] = e — 1//3 J s[q{x) , m{x)] dx , (59) 
so that we get the first order conditions 

6koP - J Sq[q,m]uk{x) dx ^ (60) 

-2(3KXkmk- Sm[q,m\uk{x)dx = Q (61) 
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We can easily check that the uniform profile q ~ q^, m ^ rfiQ (with qo and 
mo given by the a = case) is indeed a solution of these conditions. To prove 
that the uniform profile is a free energy minimum, we evaluate now the Hessian 
matrix at this point 



Ski (62) 



dqkdqi P 

-Ski (63) 



qkdrhi (3 
df 



jnrn 



drfikdrhi f3 



Ski - 2i^Afc4i . (64) 



This infinite Hessian matrix turns out to be positive definite. To prove this, it is 
sufficient to verify that each matrix Hk 

Hk = - ( - ) (65) 
+ 2K(3Xk I ^ ^ 

is positive definite, when taken at the solution point. This follows easily from 
the fact that Hq 

is positive definite (because the point considered is the solution of the a — Q 
case), and recalling that Afc < 1 for whatever k. 

We have thus proved that the local stability of the uniform profiles at a = 
are not modified in the canonical ensemble when a > 0. All canonical T(e) curves 
shown in Fig. ^ thus extend to the a > case. We now make use of the general 
discussion at the end of Section ^. We recall that results there derived state, in 
particular, that a stable canonical solution is always a stable microcanonical one, 
and that if one follows a T(e) curve, like in Fig. |[ a change of stability in the 
microcanonical ensemble may occur only where this curve has a vertical tangent 
(see [|l2| and the contribution by Chavanis |^ to this book). We thus conclude 
that the microcanonical part of Fig. ^ does not change either. This reasoning is 
however not rigorous, since we have considered only local stability; it may happen 
that when a > 0, a heterogeneous solution appears and it could yield the true 
free energy minimum or entropy maximum. Some numerical investigations we 
have carried out did not detect any such heterogeneous solution. 

We now turn to a different issue, namely the correspondence between the mi- 
croscopic models and the continuum ones. In the following we show that the Ising 
energy functional (|5^) can be derived by coarse graining from the microscopic 
Ising Hamiltonian with slowly decaying interactions 



where is the distance on a ID lattice between spins at sites i and j. The 
interaction is non integrable for a <1. The normalization N — 2"A^^^"/(1 — a) 
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ensures that the energy is extensive and the prefactor provides a convenient 
normahzation for the eigenvalues of the Hamihonian in the case of periodic 
boundary conditions. To make the correspondence between continuous and mi- 



croscopic models explicit, one has to set J = (1 — a) J/2" in Eq. 50 
Analogously, the BEG microscopic model we refer to is 

"B.O^^^Sf-^^?^. (68) 

The coarse-graining procedure we present, closely follows Rcf. |2^. We will 
keep here a heuristic level of description, although all the following calculations 
can be made rigorous using the language and techniques of large deviation the- 
ory ijl^ . A more detailed account of this derivation is given in Ref. |Q . 

The first step involves a coarse-graining procedure. We divide the lattice in 
N/n boxes, each one containing n sites. We then describe macroscopically the 
system in the Ising case by an average magnetization in each box, or, in the BEG 
case, by the average m and q. In the limit N oo, n — > oo with N/n oo, 
the system is then described by continuous functions m{x), q{x), with x G [0, 1] 
if the system length is normalized. We have to show in the following that all the 
information lost in the coarse graining procedure is unessential. 

In the mean-field case (a = 0), the Hamiltonian had a natural and exact 
expression as a function of the macroscopic parameters m (Ising) or m and q 
(BEG). The long-range Ising case has already been discussed ||2^ and it has been 
shown that the intra-box couplings can be neglected in the thermodynamic limit, 
leading to the expression for the Hamiltonian functional (js^) and for the entropy 
( ^l|) as a function of the macroscopic parameter m(x). The generalization to the 
BEG model of the method applied in Rcf. is straightforward and allows one 
to obtain the Hamiltonian functional (^7|) . The estimation is uniform over all the 
microscopic configurations p3| , and r{x, y) is the distance between the two boxes 
located in x and y. On the contrary, ii a > d, there is no way to approximate 
the Hamiltonian in the continuum limit described above as a functional of m[x) 
and q{x). 

We show now how to construct the entropy functional for the Ising model. 
Its extension to the BEG model is straightforward. We have to estimate the 
probability to obtain a certain given macroscopic configuration m{x)^ assuming 
the equiprobability of microscopic configurations. For a finite number of boxes, 
N/n, we can evaluate the probability to get a given average magnetization in 
each box 

P(mi,m2, . . . ,m„) = P{mi)P{m2) ■ ■ ■ P{mn) 

where Si is the entropy associated to the i-th box. Letting n and N/n go to 
infinity, we obtain 

P[m(a;)]~e'^r^('"(^»'^^ (70) 
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which defines the entropy functional sising of the Ising model. Let us remark 
again that large deviation techniques give a precise and rigorous meaning 
to these calculations, and allow one to extend this type of derivation to all sort 
of lattice models. The field q{x) of the BEG model can be treated in a similar 
way as m{x) and one gets the BEG entropy functional similarly. 

Chavanis (5) and, Cohen and Ispolatov |^ analyze in this book more real- 
istic off-lattice systems, self-gravitating or interacting through a 1/r" potential 
(a < 3) . They study the thermodynamic properties of these models through a 
mean-field approximation, and the phenomenology is pretty much the same as 
the one decribed here. Unfortunately, to see the analogies, one must cope with 
discrepancies in the "vocabulary" . For instance, the microcanonical first order 
transition described here is called gravitational first order transition, whereas 
the canonical first order one is called normal first order transition. Within the 
conventions of this paper one can see that, varying a control parameter (the 
interaction K for the BEG model, and a parameter controlling the short range 
cut-off in gravitational models) the system crosses over from a microcanonical 
first order transition to a microcanonical second order one (or no transition at 
all in the gravitational systems) associated with a canonical first order, and then 
recovers full ensemble equivalence, crossing the critical point. 

All these off-lattices studies raise the question of the validity of the mean- 
field approximation, which is intimately related to the scaling with TV of the 
thermodynamic variables and potentials. In our case, the mean field treatment 
is fully valid in the thermodynamic limit, but the problem was avoided using 
the physically unjustified Kac prescription Q| in which the coupling constants 
are rescaled by N. 

7 Conclusions 

In this paper a detailed comparison is made between the canonical and the mi- 
crocanonical phase diagrams of the spin-1 BEG model with mean field long range 
interactions. Since systems with long range interactions, and particularly mean 
field models, are non-additive the two ensembles need not be equivalent. The 
BEG model, for which both phase diagrams can be evaluated analytically, pro- 
vides a convenient and interesting ground for studying the distinctions between 
the two ensembles. Although the model is rather simple, it exhibits a non-trivial 
phase diagram, which in the canonical ensemble consists of both a first order 
and a second order transition lines separated by a tricritical point. It is shown 
that the microcanonical phase diagram, while it yields the same thermodynamic 
behavior in the region where the canonical transition is second order, it differs 
considerably from the canonical one in the region where the transition is first 
order. In particular it exhibits a region with a negative specific heat, a tricriti- 
cal point whose location differs from that of the canonical ensemble, and a first 
order line at which the two coexisting phases have different temperatures. The 
mechanisms which lead to these features are rather general, and thus they are 
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expected to take place in other models with long range interactions within the 
microcanonical ensemble (see also | |2^ , p7t ). 

A straightforward generalization of these results to models with algebraically 
decaying interactions are considered. It is shown that the mean field phase dia- 
gram remains valid even for these interactions, as long as they are long range. 

It would be of great interest to extend this study in the future and to an- 
alyze and compare canonical and microcanonical phase diagrams in more com- 
plex cases, where higher order multicritical points are present. Such studies may 
eventually lead to general rules which should be obeyed by microcanonical phase 
diagrams of systems with long range interactions, in analogy to the rules existing 
for systems with short range interactions, like the Gibbs phase rules, the Landau 
and Lifshitz symmetry rules for continuous transitions and others. 
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